Draft version February 2, 2008 

Preprint typeset using I4TgX style emulateapj v. 10/10/03 



CAN THE STEEP MASS PROFILE OF A1689 BE EXPLAINED BY A TRI AXIAL DARK HALO? 1 

Masamune Oguri, 2,3 Masahiro Takada, 4 Keiichi Umetsu, 5 and Tom Broadhurst 6 

Draft version February 2, 2008 

ABSTRACT 

The steep mass profile of A 1689 derived from recent detailed lensing observations is not readily reconciled 
with the low concentration halos predicted by the standard cold dark matter (CDM) model. However, halo 
triaxiality may act to bias the profile constraints derived assuming a spherically symmetric mass distribution, 
since lensing relates only to the projected mass distribution. A degree of halo triaxiality is inherent to the CDM 
structure formation, arising from the collision-less nature of the dark matter. Here we compare the CDM-based 
model predictions of triaxial halo with the precise lensing measurements of A 1689 based on the Advanced 
Camera for Surveys/Hubble Space Telescope and Subaru data, over a wide range of lOkpc < r < 2Mpc. The 
model lensing profiles cover the intrinsic spread of halo mass and shape (concentration and triaxiality), and 
are projected over all inclinations when comparing with the data. We show that the model parameters are only 
weakly constrained and strongly degenerate mainly because of the lack of information along the line of sight. In 
particular, the limits on the concentration parameter become less restrictive with increasing triaxiality. Yet, by 
comparing the obtained constraints with expected probability distributions for the axis ratio and concentration 
parameters computed from numerical simulations, we find that ~ 6% of cluster-size halos in the CDM model 
can match the A1689 lensing observations at the 2-cr level, corresponding to cases where the major-axis of 
the halo is aligned with the line of sight. Thus halo triaxiality could reduce the apparent discrepancy between 
theory and observation. This hypothesis needs to be further explored by a statistical lensing study for other 
clusters as well as by complementary three-dimensional information derived using X-ray, kinematics, and 
Sunyaev-Zel'dovich effect observations. 

Subject headings: cosmology - gravitational lensing - galaxies: clusters: individual (Abell 1689) 



1. INTRODUCTION 

The current standard model of structure formation, in which 
the universe is dominated by Cold Dark Matter (CDM), pre- 
dicts that dark matter halos have an inner cusp that is shal- 
lower than a singular isothermal sphere model. Specifically, 
the radial density profile c an be we ll fitted by an NFW pro- 
file p(r) oc r"'(l +r/r s )~ 2 ( Navarroet al.lll997h . An impor- 
tant parameter of this model is the concentration parameter 
Cyii, the ratio of the virial radius to the scale radius r s . For 
massi ve clusters, it is predicted to be ~ 4 (e.g., Bullock et al. 
1200 ll) . These are quite important predictions that should be 
confronted with observations. 

Now it is becoming possible to directly test the NFW 
predictions using gravitational lensing. A poster child ex- 
ample for this is A1689: The inner mass distribution is 
tightly constrained from more than 100 multiple images 
of background galaxies discovered from the spectacular 
deep Advanced Camera for Surveys ( ACS)ZHubble Space 
Telescope (HST) data (Broadhurs t et al.l 1200 5b, hereafter 
B05b), while the larger scale mass distribution up to the 
virial radius is obtained from the weak lensing measure- 
ments based on the wide-field S uprime-Cam/Subaru data 
llBroadhurst. Takada. Umetsu et al.ll2005al hereafter B05a). 
The two-dimensional mass profile of A 1689, reconstructed 
from the combined lensing information, does continuously 
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flatten towards the center like an NFW profile, but the fitting 
to the NFW predictions leads to a surprisingly high concentra- 
tion c v i r = 14 ± 1 .5 (B05a), compared with the expected value 
c v ii ~ 4. Such a high c oncentration is also s een in other clus- 
ters, e.g., M S2 137-23 llGavazziet al.ll2003l) and possibly CL 
0024+1654 (Kneib et al. 2003). These results may represent 
a new crisis in the current standard CDM model, offering a 
very rewarding issue to further explore. 

However, the constraints above are obtained by deproject- 
ing the two-dimensional lensing information under assump- 
tion of the spherically symmetric mass distribution. What the 
CDM model does predict is that dark halos are quite triaxial 
rather than spherical as a natural consequence of the collision- 
less dark matter nature and the fi lamentary nature of struc- 
ture f ormation (tTing & Sutoll2002l hereafter JS02; iLee et all 
2005). Therefore the results obtaine d assuming spherical 
symm etry might be biased (e.g., see iMiralda-E scude & Babull 
1995, for the s imilar discussion on A1689). In fact, 
Clowe et al. (2004) argued that the halo triaxiality affects both 
mass and concentration parameter measurements from grav- 
itational lensing using A^-body simulation results of clusters. 
Gavazzi (2005) also investigated the importance of the halo 
triaxiality and argued that a high concentration of MS2 1 37-23 
could be explained by a halo having the major axis oriented 
toward the line-of-sight. 

In this paper, we study the mass distribution of A 1689 based 
on triaxial dark halo model. Specifically, we adopt the triax- 
ial halo model of JS02, and focus on whether the apparent 
steep mass profile of A 1689 can be ascribed to triaxiality of 
the mass distribution. Throughout this paper we assume a 
concordance cosmology with the matter density Q,m = 0.3, the 
cosmological constant = 0.7, the dimensionless Hubble 
constant h = 0.7, and the normalization of matter power spec- 
trum as = 0.9. Note that one arcminute corresponds to the 
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FIG. 1. — Convergence profiles of the triaxial halo with the virial mass 
M v j, = 1O 15 A _1 M0, the concentration parameter c e = 1.15, and the triaxial 
axis ratios of a/c = 0.4 and b/c = 0.7. The halo is placed at zi = 0.3, and we 
assume the source redshift of z, s = 1-0. We consider the projection along each 
of the three principal axes: from upper to lower the dashed lines show profiles 
projected along z, y, and x (see eq. iTll. respectively. The convergence profile 
of the corresponding spherical NFW halo is also plotted by the solid line 
for comparison (see text for details). The vertical arrow indicates the virial 
radius. 



physical scale of 129/i _1 kpc for A1689 (redshift z = 0.18). 

2. A SIMPLE ESTIMATION OF THE HALO TRIAXIALITY EFFECT 
ON LENSING MEASUREMENTS 

Before going to the analysis of A 1689, we make a sim- 
ple test to demonstrate how important the halo triaxiality is 
in constraining mass profiles from a two-dimensional lensing 
m easurement. The an alysis is somewhat similar to that done 
bv lClowe et alJ ( 12004) who used high-resolution A^-body sim- 
ulations of massive clusters. Here we instead use an analytic 
model of aspherical dark halos. 

We consider a triaxial halo with the virial mass M v „ = 
lO 15 /z _1 M0 7 , placed at z\ = 0.3, and adopt the model mass 
profile given in JS02: 
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We adopt typical model parameters for a halo of \Q l5 h~ l MQ\ 
the triaxial axis ratios are a/c = 0.4 and b/c = 0.7, and 
the concentration parameter c e = R e /Ro, where R e is de- 
fined such that the mean density enclosed within the ellip- 
soid of the major axis radius R e is A e il(z)p a -it(z) with A e = 

5 A v i r (c 2 /ab) , is chosen to be c e = 1.15. We have checked 
that the spherically-averaged radial mass profile of the triax- 
ial halo is quite similar to the spherical NFW profile that is 
specified by the virial radius r v i r =R e /QA5, as proposed in 
JS02 (see Figure 14 in JS02), and the concentration parame- 
ter c v [ r = 4 s . However, it is non-trivial for these triaxial and 

7 The virial mass is defined by spherically averaging the halo mass dis- 
tribution (the triaxial mass profile for our case) around the halo center and 
then by finding the sphere inside which the mean overdensity reaches A v ; r 
predicted in the top-hat spherical collapse model. 

8 Since the overdensity in the triaxial model, A c , is at least 5 times larger 




FIG. 2. — Constraint contours in the virial mass and halo concentration 
parameter space, obtained by fitting the mock data of triaxial halos to the 
spherical NFW halo model. The contours show 68%, 95%, 99.7% confidence 
limits (corresponding to A\ 2 = 2.3, 6. 17 and 11.8, respectively). From left to 
right, the constraint contours from the convergence profiles of the triaxial halo 
projected along the principal axes x, y, and z (as in FigureQ, respectively, 
are shown. For comparison, the square symbol shows the best-fitting model 
for the convergence profile obtained by projecting the a priori spherically- 
averaged mass profile of the input triaxial halo. 



spherical models whether to yield similar lensing maps as a 
result of the line-of-sight projection 9 . To make this clear, Fig- 
ure \l\ compares the circularly-averaged convergence profiles 
for the spherical and triaxial halos. For the triaxial halo, we 
consider the projection along each of the three principal axes. 
It is clear that the surface mass density of the triaxial halo 
depends strongly on the projection direction. Therefore it is 
quite likely that adopting a spherical halo model causes a bias 
in estimating the mass and profile parameters for an individual 
cluster in reality. 

To see this more clearly, we perform the following test. 
First, we generate an "observed" surface mass density pro- 
file: We consider 20 bins logarithmically spacing over the 
range r = [10~ 2 , l]/i _1 Mpc, and generate the convergence pro- 
file At(r), where the mean value for each bin is taken from the 
triaxial halo model and the Gaussian random error of stan- 
dard deviation A(log 10 K) = 0.1 is added to each bin. Then, 
assuming the spherical NFW density profile, we constrain the 
virial mass (M v ; r ) and halo concentration parameter (c v ir) by 
fitting the model predictions to the "observed" profile. The 
constraint contours in the M wu -c v ; r plane are shown in Figure 
13 demonstrating that the best-fit parameters depend strongly 
on the projection direction. For example, the convergence 
profile projected along the major (minor) axis yields a sig- 
nificant overestimation (underestimation) by 20-30% in both 
the mass and concentration parameters. It should be noted 
that the bias direction is orthogonal to the degeneracy direc- 
tion of the error ellipse, implying the systematics is very im- 

than the spherical overdensity A v i, , the concentration parameter in the triaxial 
model tends to be smaller than those in the spherical model (approximately 
given as c e f» 0.45c v j r ). See JS02 for details. 

9 The lensing convergence field re(r) is given in terms of the surface mass 
density S(r) as n(v) = E(r)/£ cl , where S cl is the lensing critical density 
specified for a background cosmology and lens and source redshifts (see 
Schneider e t all 19921) . 
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Table 1. Best-fit parameters for the model fitting to A 1689 



X 2 R c [Mpc] 


C e 


a/c 


b/c 


M vir [lO 15 M ] 


Xobs LSO -0.55 
v 2 +v 2 i 5(Y*0.11 


7 6 +L0 


0.1« 0.75) 

45+0.07 


0.85(> 0.32) 
6+ - 20 


1 8+0-4 

1 7 +0 ' 3 
l./_0 2 



NOTE. — Best-fit parameters and the 1-c errors are presented, for the triaxial halo model fitting to the A1689 lensing information. Note that the virial mass 
Mvjj is not a model parameter, but is derived from the constrained model parameters R e , c e , a/c, and b/c (see text for details). The lower row shows the results 
when the theoretically expected probability distributions for the axis ratios are employed as the priors of the fitting. For some parameters, only a lower or upper 
limit is obtained. 



portant. In fact, we fail to recover the model parameters of 
the spherically-averaged triaxial profile at more than 3-a level 
when the halo is projected along the major- or minor-axis di- 
rections. 

3. APPLICATION TO A 1689 

In this section, we indeed make a quantitative interpretation 
of the lensing measurements of A 1689 in terms of the triax- 
ial halo model in order to derive more reliable constraints on 
the three-dimensional mass distribution than does assuming 
the spherical halo model. Our particular attention is paid to 
whether or not the steep mass profile claimed in B05a can be 
reconciled with the CDM-bas ed triaxial h a lo mo del of JS02. 

Following the method of Osuri et al. (2003), the model 
lensing convergence field for a triaxial halo is specified by 
7 parameters. As explained around equation Q, the virial 
radius of the triaxial halo modeling R e , the concentration pa- 
rameter c e , and the axis ratio parameters a/c and b/c are con- 
sidered. Further, we adopt three angle para meters to specify 
the halo orientation: 9, qt> (see Figure 1 in lQguri et alJl2003l) 
and the third angle between major axis of an ellipse of the 
projected mass distribution and the RA direction on the sky. 
Note that the parameter 9 gives the angle between the major 
axis of the triaxial halo and the line-of-sight direction. We 
again note that, in this model, the virial mass M v „ is not a free 
parameter but is derivable from the model parameters (R e , c e , 
a/c, b/c). Throughout this paper, we assume (z s ) = 1 for the 
mean redshift of source galaxies as done in B05a. 

We then constrain the 7 model parameters from x 2 fitting to 
the observed convergence map obtained from the ACS/HST 
and Suprime-Cam/Subaru data (e.g. see B05a for the similar 
method): 

Xobs = XHST + XSubara ' O) 

where we adopt flat priors of 0.1 < a/c(< 1) and 0.1 < 
b/c(< 1) for the axis ratios because a small axis ratio such 
as a/c < 0.1 is unrealistic for a virialized halo. For the ACS 
data, we use the circularly-averaged profile of the obtained 
convergence map that is given in 12 bins linearly spacing over 
the range r = [0f038,0.'97] (see Figure 22 in B05b or Figure 3 
in B05a) 10 . Hence, the Xhst lS given by 



Xhst _ 
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where k™ is the model prediction of the triaxial halo for the 
z'-th radial bin, computed as in Figure and the k, and er,- are 

10 It is not straightforward to use the two-dimensional mass map obtained 
from the ACS strong lensing analysis because of the complex, non-linear 
error propagation. The bin width of the one-dimensional convergence profile 
we use here is sufficiently broad to ensure that the errors between different 
bins are independent (B05b). 



the observed value and 1-er error, respectively. As for the Sub- 
aru data, we use the two-dimensional convergence map k(t) 
reconstructed from the joint measurements of the weak lens- 
ing distortion and magnification bias effects on background 
galaxies (Umetsu et al. 2005, in preparation). The map is 
given on 21 x 17(= 357) grids for the area of w 30' x 24'. 
Note that the wide-field Subaru data allows us to probe the 
mass distribution on larger scales > 1' covering up to the 
virial radius ~ 20' (B05a), thus ensuring that the ACS and 
Subaru information are independent. The pixel-pixel co- 
variance matrix defined as Vy = (<5K(r,) c>K(r,)} can be ac- 
curately estimated based on the maximum likelihood map- 
making method, assuming that the error of weak lensing dis- 
tortion arises from the random intrinsic ellipticities of back- 
ground galaxies, while the magnificatio n bias error is th e Pois- 
sonian noise of the number counts ( Schneider et al. 2000; see 
Umetsu et al. in preparation for the details). The x| ubaru can 
be therefore expressed as 

(5) 



xLaru = E [« m (r«)-«(r,)] (V-% [ K m (rj)- K (rj) 

i-j 



where V -1 is the inverse of the covariance matrix. 

Tabled lists the best-fit values for the 4 model parameters 
of our interest as well as the virial mass derived, with the l-a 
error including the other 6 parameter uncertainties (estimated 
from A% 2 = x 2 ~ Xmin = 1 m tne 7 parameter space). The 
minimum \ 2 is Xmin/dof = 378/362 for 362 degrees of free- 
dom, thus showing an acceptable fit. It is clear that the halo 
triaxiality significantly weakens constraints on the halo con- 
centration parameter c e = 7.6^ g, compared to the constraint 
c v j r = 14 ± 1 .5 (c e w 0.45c v j r ) derived when a priori assuming 
the spherical NFW model (B05a). In fact, a smaller concen- 
tration compatible with the theoretical expectation c e = 1-2 
is marginally allowed at 1-cr level. 

This is more clearly explored in Figure |3j which shows 
the error contours in the two parameter subspace. We also 
show the error contours in the c e -M wu plane in Figure |4] to 
compare our result with Figure 3 of B05a. Apparently, the 
parameter constraints are significantly degenerate. For ex- 
ample, as a triaxial halo with smaller a/c is considered, a 
broader range of the halo concentration is allowed, imply- 
ing that the lensing strength is as quite sensitive to the halo 
triaxiality as is to the halo concentration c e . However, the 
axis ratios can be hardly constrained, obviously because of 
a lack of the halo shape information along the line-of-sight: 
Any values of a/c and b/c are allowed to within 3-a level, 
though profiles close to spherical models (a/c ~ 1) are rela- 
tively disfavored because the Subaru convergence map is not 
circularly symmetric. More intriguing is that the degeneracy 
direction in the c e -M vlr plane (Figure|4j is almost orthogonal 
to that for the spherical case (Figure 3 in B05a), as illustrated 
in Figure |2 As a result, the mass constraint can be recon- 
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FIG. 3. — Projected constraint contours in the two parameter subspace from the 7 parameter space for the triaxial halo model, matching the two-dimensional 
mass distribution of A1689 reconstructed from the ACS/HST and Subaru data. The dark shaded regions show the 68%, 95% and 99.7% confidence intervals 
corresponding to A\ 2 = 2.3, 6.17, and 11.8 in the 7 parameter space, respectively. It is clear that the parameter constraints are significantly weakened by 
inclusion of the halo triaxiality and show strong degeneracies in each parameter space. Note that the concentration parameter of the triaxial halo model, c e , is 
approximately related to that of the spherical model via c e ~ 0.45c v j r (JS02). The solid contours show the improved constraints when the CDM-based theoretical 
predictions of probability distributions for the axis ratios a/c and b/c are added as the prior to the x 2 -ntting (see text for details). However, a broad range of the 
parameters are still allowed mainly because of a lack of the halo shape information along the line-of-sight. 



ciled with the mass estimate derived from th e X-ray ob serva- 
tion (w 1.0 x 10 15 M Q in lAndersson & Madeiskl2 004) at 3-er 
level. Thus the halo triaxiality co uld resolve the mass dis- 
crepancy reported in the literature ( Miralda-Escude & Babul 
119951 lAndersson & Madelskll2004 B05b), as we will again 
discuss later in more detail. 

To improve the weak constraints above, we employ the 
prior knowledge of what kinds of halo shapes are expected 
for cluster-scale halos within the CDM clustering scenario. 
We use the probability distribution functions (PDFs) of the 
axis ratios, which are derived in JS02 using the A/-body simu- 



lations. The PDFs are given as a function of redshift and halo 
mass; a cluster-scale halo at z ~ is typically characterized 
by a triaxial shape with the mean value a/c ks 0.45 and the 
standard deviation a(a/c) ~ 0.1. For simplicity, we add the 
following Gaussian prior to the \ 2 ( e q- G)) : 

Xl = -2.01n [p(a/c)p(a/b\a/c)] , (6) 

where the subscript 'ar' stands for axis ratios and the PDFs 
p(a/c) and p(a/b\a/c) are given by equations (17) and (19) 
in JS02, respectively. The resulting constraints are shown in 
Table [2 and Figures [3] and 0] However, the constraints are 
only weakly tightened and the main results we have so far 
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FIG. 4. — As in the previous figure, the constraint contours are shown in 
the c e - M v j, plane for the triaxial halo model, where the virial mass M viI is 
derived from the model parameters constrained in Figurel3l Note that the con- 
centration parameter of the triaxial halo model, c e , is approximately related 
to that of the spherical model via c e ~ 0.45c v j r (JS02). 
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FIG. 5. — Constraint contours (solid lines) in the c e -a/c plane as in the 
previous figure, but the radius R e is fixed to R e = 1 .4Mpc which is consistent 
with observational constraints at 1 — a level (see Table [TV Note that these 
results include the uncertainties of the other parameters besides R e . The con- 
straints are directly compared with the theoretical predictions (dashed lines) 
giving the 68%, 95% and 99.7% probabilities for the c e and a/c distribution 
derived using the model of JS02. While the PDF of c e slightly depends on 
halo mass, in computing the theoretical predictions we fix the virial mass to 
M V11 = 1.5 X 1O 15 M which roughly corresponds to R e = 1 .4Mpc for a given 
range of (c e , a/c, b/c) we have considered here. It is clear that the observation 
and theory are consistent at ~ 2-a level (see text for more details). 



found are not largely changed: a broad range of the halo con- 
centration or the virial mass is still allowed. 

Hence, to more highlight the constraint on the halo con- 
centration, Figure [5] shows the error contours in the c e -a/c 
plane, when the radius R e is fixed to R e = 1 .4Mpc which is 
consistent with observational constraints at 1-er level. Note 
that the other parameters are allowed to vary. The clear trend, 
a smaller c e is favored with decreasing a /c, is apparent. This 
result allows a more direct comparison with the theoretical 




cos 9 

FIG. 6. — The likelihood function against the angle parameter 6 between the 
major axis of triaxial halo and the line-of-sight direction, assuming the priors 
for the PDFs of axis ratios and halo concentration. The horizontal dotted lines 
show 68%, 95%, 99.7% confidence limits (Ax 2 = 1, 4 and 9, respectively). 
It is clear that an alignment of the halo elongation along the line-of-sight is 
favored in order for the CDM-based halo model to match the A 1689 lensing 
measurements. 



prediction, since JS02 suggested that the halo concentration 
c e be tightly correlated with the axis ratio a/c. The PDF of c e , 
p(c P ), is given by equation (20) in JS02 (also see Cteuri et al. 
2003). The PDF is approximated with a log-normal distribu- 
tion, where the mean value of c e depends on halo mass, a/c 
and redshift (see eq. [23] in JS02) and the standard deviation 
<r(\Tic e ) « 0.3. Combining with the PDFs (eq. @) for the axis 
ratios, the dashed contours in Figure[5]show the 1-, 2- and 3- 
a probabilities for the c e and a/c parameters, expected for 
a halo at z = 0.18, which are estimated from + Xce = 2-3, 
6.17, and 1 1 .8, respectively, where xl e = _ 2 In [p(c e )] (the sub- 
script 'ce' stands for the concentration parameter c e ). It is 
clear that the observation and theory overlap at ~ 2-a level. 
To see this more quantitatively, we compute the fraction of ha- 
los whose axis ratio a /c and concentration parameter c e fall 
inside the contours constrained from the observation. We find 
that, if we take 1-, 2-, and 3-er limits as the observational con- 
straints, then the fraction becomes 2.7%, 6.1%, and 11.6%, 
respectively. In contrast, if we a priori assume the spheri- 
cal halo model (B05a), the fraction is significantly reduced to 
< 0.1%, even if we take the 3-cr limit as observational con- 
straints. Hence, we conclude that the A1689 lensing mea- 
surements are indeed compatible with the CDM-based triaxial 
halo model, if the A1689 represents a rare population (« 6% 
fraction) of cluster-scale halos (also see the discussion in the 
next section). 

Finally, we look at the constraints on the angle parameter 
8 between the major axis and the line-of-sight direction. Fig- 
ure [6] shows the likelihood function against 8, derived from 
the observational constraints x^ bs plus the theoretical PDFs of 
axis ratios and concentration parameter, xlt + xl e - O ne can 
find cos6> > 0.8 (8 < 37°) is obtained at 2-a level. Thus the 
major axis of the triaxial halo, when the halo has a reasonable 
concentration of c e = 1-2, is favored to be aligned with the 
line-of-sight direction in order to match the A 1689 steep mass 
profile. 

4. SUMMARY AND DISCUSSION 
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In this paper, we have employed the triaxial halo model 
of JS02 to extract more reliable information on the three- 
dimensional mass distribution of the massive cluster A 1689, 
from the precise lensing measurements done based on the 
ACS and Subaru data. We have shown several important re- 
sults. First, based on the simple thought experiment in ^2] it 
was demonstrated that the halo triaxiality could cause a sig- 
nificant bias in estimating the virial mass and concentration 
parameter from the lensing information, if a spherical halo 
model is a priori assumed for the model fitting (see Figure 
|2jl. In particular, both the virial mass and concentration pa- 
rameter are overestimated (underestimated) when the triax- 
ial halo is projected along the major (minor) axis. Second, 
we derived constraints on the parameters of the triaxial halo 
model from the fitting to the lensing measurements of A 1689, 
by changing all 7 parameters that specify the triaxial model; 
the virial radius and concentration parameter of the triaxial 
modeling, two axis ratio parameters, and three angle parame- 
ters that determine the halo orientation. We have shown that 
the halo triaxiality significantly weakens the parameter con- 
straints, compared to those derived when a priori assuming 
the spherical halo model. In addition, the parameter con- 
straints are strongly degenerate, mainly because of a lack of 
the halo shape information along the line-of-sight (see Table 
[2 and Figures |3] and |4}. As a result, the triaxial halo model 
with a theoretically expected halo concentration c e = 1 — 2 is 
acceptable at l-er level, while the halo shape parameters (axis 
ratios and orientation) are hardly constrained. To further im- 
prove the model constraints, we used the theoretical predic- 
tions of probability distributions for the halo shapes, derived 
in JS02, as the priors to the model fitting (see Figures [3] H E] 
and|6j. Most interestingly, we found that about 6% population 
of cluster-scale halos, expected from the CDM structure for- 
mation, can match the A 1689 mass profile, if the major axis 
is almost aligned with the line-of-sight direction. Hence, al- 
though it has been shown that the triaxiality has a great impact 
on strong lensing probabiliti es (Oguri et al. 2003: [Dalai et alJ 
2004: iQguri & Keetonl2 004). our results clearly demonstrate 
that the halo triaxiality is also very important for exploring 
the three-dimensional mass distribution of an individual clus- 
ter from the lensing measurements. 

It is encouraging that the triaxial halo model yields the 
lensing-constrained mass of A 1689 consistent with the mass 
estimate derived from the X-ray observation at 3-a level 
(see Figure |4j, while the significa nt discrepancy of factor 2 
has been reported in the literature ( Miralda-Escu de & Babull 
H995U Andersson & Madeiski 2004; B05b). We have checked 
that the consistent mass estimate appears when the major axis 
of the triaxial halo is oriented along the line-of-sight direc- 
tion, as inferred from Figure |2] Note that the X-ray mass es- 
timate is likely to be less af fected by the ha lo triaxiality than 
the lensing, as described in Gavazzi (2005). Therefore, it is 
interesting to notice that the triaxial halo elongated along the 
line-of-sight could resolve both the previously reported dis- 
crepancies of the high concentration as well as of the lensing 
and X-ray mass estimates. Unfortunately, however, this hy- 
pothesis cannot be further tested by the lensing information 
alone we have used, because the halo shape parameters (axis 
ratios and orientation) could not be well constrained. Improv- 
ing the halo shape constraints will allow us to make a more 
severe, quantitative test of the triaxial halo model, as it is clear 



from Figure|5] This would be possible by combining the lens- 
ing measurements with the detailed observations of kinemat- 
ics, X-ray emission and/or the Sunyaev-Zel'dovich (SZ) ef- 
fect, since projections affect the observables in different ways 
(e.g., iGirardi et alJll997l: IZaroubi et alJl2001t iMarshal et alJ 
2003). For instance, as noted above, the mass profile re- 
constructed from the X-ray surface brightness profile is much 
less sensitive to the halo orientation, thus by comparing two- 
dimensional profiles from X-ray and lensing we will be able 
to obtain information on the degree of elongation along line- 
of-sight. The SZ effect is also quite usefu l since its projectio n 
effect is complementary to that of X-ray (Lee & Suto 2004). 

While our triaxial modeling is much more realistic than the 
simple spherical fit, there is still room for improvement. A 
caveat is that the PDFs of axis ratios and concentration (JS02) 
were derived for less massive clusters (M vir < 1O 14 M ). Thus 
it is unclear if the extrapolated predictions still hold for very 
massive clusters (M vn - > 10 15 M Q ) considered here. More im- 
portantly, the selection effect, which we have neglected in 
this paper, might be important. A 1689 is well known as 
the strong-lensing cluster, with the largest Einstein radius (~ 
50")- Strong lensing probabilities are very sensitive to shapes 
and orientations of dark halos such that a halo with large triax- 
iality, high concentration and the major axis aligned with line- 
of-sight has higher lensing pr obabilities (Oguri et alJl2003t 
Dalai et alJl2004t IQ guri & Keetoi] 120041) . Therefore A1689 
might be indeed such a cluster, as implied in this paper. Thus 
taking into account the selectio n effect could further reduce 
the discrepancy of A1689 (e.g. JHennawi et aTll2005). How- 
ever, in practice it is quite difficult to quantify the selection ef- 
fect. For strong lens systems discovered in a well-defined ho- 
mogeneous survey, s uch as the large-separation lensed qua sar 
SDSS J1004+4112 foiada et alJl2003l: TOguri et ai1 l200l . it 
will be possible to take this selection effect into account in 
a correct manner. 

Another important implication we have found is that the 
halo triaxiality can be a source of systematic errors in esti- 
mating the virial mass from the lensing measurement. While 
it has been shown that the lensing-based mass determina- 
tion suffers from the projection effect due to an interven- 
ing m atter, not necessarily associated with a clust er of in- 
terest (Metzleret al. 2001; Wambsgans s et alJl2005t) . our re- 
sults imply that the hal o shape itself causes a bias in the mass 
determination (also see Bartelman nl 1 9951 IClowe et al.ll2004t 
IHamana et al l200l . The accurate mass determination is cru- 
cial to achieve the full potential of future cluster cosmology 
experiments such as expected from Planck, where the clus- 
ter redshift distribu tion is measu red to constrain cosmologi- 
cal models (Haimanet al. 2001). Therefore, it is again very 
interesting to address how the mass determination can be im- 
proved for an individual cluster as well as in a statistical sense 
by combining the lensing, X-ray and SZ effect measurements. 
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